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SECTION  I 


INTRODUCTION 


Engineering  techniques  for  predicting  and  measuring  stability  and  controlla- 
bility of  flexible  airplanes  are  in  a state  of  rapid  development.  Increased 
attention  to  the  problem  results  from  complexity  of  current  design  require- 
ments, promise  of  future  performance  improvements,  and  recent  theoretical 
advances . 

Historically,  interaction  of  structural  flexibility  with  flight  modes  has 
been  important  for  very  large  airplanes,  such  as  the  B-52.  New  design 
requirements,  especially  flight  at  high  dynamic  pressures,  are  increasing  the 
severity  of  the  effects  so  that  smaller  aircraft  now  must  be  aeroelastically 
designed  and  tested. 

Modern  automatic  stability  aid  control  augmentation  systems  are  incorporated 
on  flexible  airplanes  only  after  thorough  design  and  testing,  often  with 
testing  revealing  design  deficiencies.  The  utility  of  more  complex  systems  is 
readily  foreseen.  A promising  use  of  c ntrol  systems  is  to  provide  stability 
to  the  extent  that  the  configurator  may  select  a low  drag  or  highly  maneuverable 
design  that  would  otherwise  be  unflyable.  The  same  systems  can  also  control 
elastic  modes  to  relieve  stresses  and  prevent  flutter.  This  is  the  Control 
Configured  Vehicle  (CCV)  concept.  Substantial  cost  savings  and  performance 
Improvements  can  be  predicted  for  this  approach.  However,  the  economic 
risk  of  constructing  such  a vehicle  can  be  justified  only  if  the  analysis 
methods  used  to  predict  aeroelastic  controllability  have  been  proven  accurate. 

A theoretical  development  which  has  provided  the  analytical  accuracy  needed 
is  the  finite  element  approach  to  numerical  solution  of  linear  aerodynamic 
theory,  both  steady  and  unsteady.  An  example  of  a computer  program  system 
which  incorporates  this  type  of  analytical  aerodynamics  is  Level  2 FLEXSTAB 
(Ref.  1).  This  program  will  require  verification  using  actual  flight  test 
data  of  aircraft  with  significant  control  and  elastic  mode  interaction.  The 
specific  type  of  flight  testing  required  and  the  method  of  verification  are 
themselves  subjects  of  theoretical  interest. 

A comprehensive  theory  of  testing,  known  generally  as  "parameter  estimation," 
has  been  developing.  Reference  2 is  a compendium  of  current  approaches  to 
extending  this  theory  to  aircraft  testing.  The  procedure  is  to  record 
aircraft  control  inputs  and  responses  during  maneuvers,  and  then  use  the 
measured  data  to  compute  refinements  to  a theoretical  model  of  the  aircraft. 
Application  of  these  methods  to  flexible  aircraft,  including  coupling  of 
flight  and  vibration  modes,  has  not  been  attempted  because  of  a lack  of 
suitable  quality  flight  test  data. 

A suitable  set  of  flight  test  data  for  verification  of  predictive  analyses  and 
development  of  parameter  estimation  techniques  now  exists,  created  during  B-52 
CCV  program  flight  testing  completed  in  December  1973.  Analytical  data  required 
for  parameter  estimation  and  for  verification  of  program  systems  such  as 
FLEXSTAB  2 are  also  available.  Since  the  tests  were  not  conducted  for  these 
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specific  purposes,  a considerable  amount  of  data  processing  and  qualification 
has  been  done. 


1 . 1 Background 

The  Boeing  Company,  Wichita  Division,  has  conducted,  in  conjunction 
with  AFFDL,  under  Contract  F33615-71-C-1926,  the  B-52  phase  of  the  CCV 
program.  This  program  was  begun  in  1971,  flight  testing  was  completed 
by  the  end  of  1973,  and  final  reports  were  released  early  in  1975. 

The  program  included  flight  demonstration  of  five  CCV  concepts  — 

Ride  Control,  Maneuver  Load  Control,  Gust-induced  Structural  fatigue 
Relief,  Flutter  Mode  Control,  and  Augmented  Stability. 

Although  the  use  of  the  data  for  parameter  estimation  was  not  foreseen, 
the  large  number  of  stress  and  motion  sensors  used  and  recorded,  and 
the  type  of  maneuvers  performed  to  demonstrate  the  systems  provided  a 
set  of  data  highly  suitable  for  parameter  estimation  study.  The  data 
existed  on  the  raw  flight  test  tapes,  't  had  not  heretofore  been  put 
in  a suitable  form.  Data  from  most  of  the  sensors  were  processed 
during  the  CCV  program  for  low  level  turbulence  flights  — but  only  a 
few  items  were  reduced  and  documented  for  the  maneuvering  flights. 

In  addition,  theoretical  models  of  the  B-52E  using  state-of-the-art 
aerodynamic  numerical  methods  were  developed  during  the  CCV  program. 
These  were  of  the  type  required  for  parameter  estimation  studies. 

Data  of  the  form  required  to  verify  predictive  programs  such  as 
FLEXSTAB  2 were  also  available. 

References  3,  4,  and  5 document  the  parts  of  the  previous  work 
applicable  to  this  program. 

1.2  Program  Objective 


This  program  provides  a set  of  B-52E  test  data,  analytical  model  of 
the  test  airplane,  and  a preliminary  parameter  estimation  analysis  of 
the  data  using  a method  suitable  for  multiple  mode  elastic  airplanes. 


The  test  data  are  on  digital  computer  tapes  in  a form  suitable  for 
further  parameter  estimation  studies.  Airplane  descriptive  data  are 
sufficient  so  that  the  flight  maneuvers  can  be  modelled  with  programs 
such  as  FLEXSTAB  2.  Test  and  analytical  responses  are  compared  as 
part  of  the  qualification  process.  The  cross- correlation  parameter 
estimation  has  extracted  eigenvalues  for  the  flight  modes  and  several 
flexible  modes  from  the  test  data.  Recommendations  are  made  concerning 
use  of  this  analysis  method. 
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SUMMARY 


Flight  test  data  and  a test  vehicle  description  are  presented  for  dual 
purposes;  (a)  to  provide  a means  of  verifying  predictive  flexible  airplane 
response  computer  programs,  and  (b)  to  provide  data  for  developing  methods  of 
parameter  estimation  suitable  for  such  aircraft.  The  tests  were  conducted 
during  the  AFFDL-Boeing  CCV  program  using  the  Air  Force  Flight  Research 
airplane  NB52E  AF56-632. 

The  test  data  for  the  five  transient  response  maneuvers  and  three  control 
surface  frequency  sweep  conditions  are  available  on  digital  computer  tapes. 
Comparisons  with  theoretical  responses  are  made  to  verify  both  the  airplane 
descriptive  data  and  the  test  data  calibration  and  digitizing  process. 

Eigenvalues  for  rigid-body  and  six  elastic  modes  were  obtained  from  symmetric 
aileron  sweep  test  data,  showing  excellent  agreement  with  theoretical 
eigenvalues.  Transient  responses  to  aileron  motion  computed  using  equation 
coefficients  derived  from  the  sweep  test  data  compared  favorable  with  actual 
responses.  Consistent  eigenvalues  were  not  obtained  from  the  transient  test 
data. 

The  parameter  estimation  conducted  introduces  techniques  for  selecting  and 
measuring  state  vectors,  for  permitting  variation  of  rate  and  displacement 
aerodynamic  derivatives  with  frequency,  and  for  a specialized  conversion  of 
difference  equation  coefficients  to  differential  equation  coefficients.  These 
techniques  are  recommended  to  future  investigators  of  parameter  estimation 
for  flexible  airplanes.  In  addition,  the  value  of  structural  deformation 
Instrumentation  (i.e.,  strain  gages)  is  emphasized. 
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SECTION  III 

TEST  VEHICLE  CONFIGURATION 


The  teet  vehicle  for  the  Control  Configured  Vehicles  program  was  an  Air  Force 
flight  research  airplane,  NB-52E  AF5tf-632.  This  airplane  was  previously  used  as 
the  test  vehicle  for  the  Load  Alleviation  and  Mode  Stabilization  (LAMS)  program. 

This  airplane  was  chosen  for  the  CCV  program  because  it  represents  a class  of 
large  flexible  aircraft  which  could  be  expected  to  benefit  greatly  from  the 
concepts  to  be  demonstrated.  The  airplane  has  a wing  span  of  185  feet,  fuselage 
length  of  157  feet  and  fin  height  of  48  feet.  Maximum  gross  weight  is  approxi- 
mately 450,000  pounds  and  for  this  program  a range  of  weights  from  242,000  to 
373,000  pounds  was  utilized, 

The  data  required  for  the  application  of  parameter  estimation  methods  were 
selected  from  the  light  grous  weight  configuration,  The  gross  weights  for 
the  eight  flight  tests  selected  ranged  from  242,000  m 272,300  pounds. 

3.1  Test  Vehicle  Modification 

Modifications  previously  incorporated  for  the  LAMS  program  included  instal- 
lation of  two  Electronic  Associates,  Inc.  Model  TR-48  analog  computers  in 
the  lower  nose  section,  installation  of  integrated  electrohydraulic  servo- 
actuators  on  existing  elevators,  rudder,  inboard  ailerons  and  spoilers, 
conversion  of  the  pilot's  control  system  to  fully  Fly-By-Wire  (FBW) , and 
installation  of  extensive  interface  electronics  and  instrumentation. 

The  major  CCV  modifications  include  the  addition  of  two  horizontal  and  one 
vertical  canard  for  the  ride  control  system,  new  three-Begment  flaperons  and 
an  outboard  aileron  for  the  maneuver  load  and  flutter  mode  control  systems, 
and  2,000  pounds  of  lead  ballast  in  the  nose  of  each  external  fuel  tank  to 
create  a wing  flutter  mode  within  the  operational  envelope  of  the  airplane. 
The  new  and  existing  surfaces  and  ballast  locations  are  shown  in  Figure  1. 

A detailed  discussion  of  airplane  modifications  is  presented  in 
Reference  6. 

3.2  Instrumentation 

The  measurement  station  locations  of  the  data  acquired  for  this  program 
are  shown  in  Figures  2 and  3.  The  locations  of  the  gust  probe,  acceler- 
ometer and  rate  gyros  are  shown  in  Figure  2.  Bending  moments  were 
recorded  at  the  locations  shown  in  Figure  3, 

In  addition  to  the  measurements  shown  in  Figures  2 and  3,  there  are 
position  indicator  recordings  for  the  control  surfaces  shown  in  Figure  1. 

A complete  list  of  recorded  measurements  is  shown  in  Table  1.  The 
measurements  are  shown  in  the  order  used  in  transmitting  the  flight  test 
data  on  the  magnetic  data  tapes.  The  data  tapes  are  described  in 
Appendix  C. 
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Figure  1.  CCV  Flight  Control  Surfaces 
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Figure  2.  Sensor  Locations 


TABLE  1.  RECORDED  FLIGHT  TEST  MEASUREMENTS 


NO. 


ITEM 


LOCATION  UNITS 


I 

Vertical  Bending  Moment 

BS760 

2 

Vertical  Bending  Moment 

BS1222 

3 

Vertical  Bending  Moment 

BS1412 

4 

Lateral  Bending  Moment 

BS760 

5 

Lateral  Bending  Moment 

BS1222 

6 

Lateral  Bending  Moment 

BS1412 

7 

Vertical  Bending  Moment 

LWS222 

8 

Vertical  Bending  Moment 

LWS820 

9 

Vertical  Bending  Moment 

LWS974 

10 

Chordwise  Bending  Moment 

LWS222 

11 

Chordwise  Bending  Moment 

LWS820 

12 

Chordwise  Bending  Moment 

LWS974 

13 

Vertical  Bending  Moment 

RWS222 

14 

Vertical  Bending  Moment 

RWS820 

15 

Vertical  Bending  Moment 

RWS974 

16 

Chordwise  Bending  Moment 

RWS222 

17 

Chordwise  Bending  Moment 

RWS820 

18 

Chordwise  Bending  Moment 

RW3974 

19 

Vertical  Bending  Moment 

LHTBL56 

20 

Vertical  Bending  Moment 

RHTBL56 

21 

Vertical  Bending  Moment 

FS135 

22 

Pitch  Rate 

BS860 

23 

Roll  Rate 

BS860 

24 

Yaw  Rate 

BS860 

25 

Vertical  Acceleration 

BS172 

26 

Vertical  Acceleration 

BS860 

27 

Vertical  Acceleration 

BS1655 

28 

Lateral  Acceleration 

BS172 

29 

Lateral  Acceleration 

BS860 

30 

Lateral  Acceleration 

BS1655 

31 

Vertical  Acceleration 

LWBL565 

32 

Vertical  Acceleration 

LWBL925 

33 

tr  - — a.  J A amaI  a**«**4  nn 

VCLLlLai 

LWS1359 

34 

Vertical  Acceleration 

Left  Ext.  Tank,  Nose 

35 

Vertical  Acceleration 

RWBL565 

36 

Vertical  Acceleration 

RWBL925 

37 

Vertical  Acceleration 

RWS1359 

38 

Vertical  Acceleration 

Right  Ext.  Tank,  Nose 

39 

Vertical  Acceleration 

Left  Inboard  Nacelle 

40 

Lateral  Acceleration 

Left  Inboard  Nacelle 

41 

Vertical  Acceleration 

Left  Outboard  Nacelle 

42 

Lateral  Acceleration 

Left  Outboard  Nacelle 

43 

Angular  Acceleration,  Y-Axis 

LWS540 

44 

Angular  Acceleration,  Y-Axis 

LWS1359 

45 

Angular  Acceleration,  Y-Axis 

RWSJ359 

In. -LBS 


In. -LBS 
Deg/Sec 
Deg/Sec 
Deg/Sec 


Rad/Sec^ 

Rad/Sec2 

Rad/Sec^ 
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RECORDED  FLIGHT  TEST  MEASUREMENTS  (Concluded) 
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NO. 

ITEM 

LOCATION 

UNITS 

46 

Delta  Inboard  Aileron 



Deg 

47 

Delta  Rudder 

- 

k 

48 

Delta  Elevator 

- 

49 

Delta  Vertical  Canard 

- 

50 

Delta  Left  Horizontal  Canard 

- 

51 

Delta  Right  Horizontal  Canard 

- 

52 

Delta  Outboard  Aileron 

- 

53 

Delta  Outboard  Flaperon 

- 

54 

Delta  Inboard  Flaperon 

- 

55 

Angle  of  Attack 

Gust  Boom 

r 

56 

Angle  of  Side  Slip 

Gust  Boom 

Df 

ig 

57 

Forward  Velocity 

Gust  Boom 

Ft/Sec 

MEASUREMENT  SIGN  CONVENTION 


Vertical  Bending  Moments 

Lateral  Bending  Moment 

Chordwlse  Bending  Moment 
Pitch  Rate 
Roll  Rate 
Yaw  Rate 

Vertical  Acceleration 
Lateral  Acceleration 
Angular  Acceleration,  Y-Axis 
Control  Surface  Deflection 
Angle  of  Attack 
Angle  of  Side  Slip 


Fuselage,  + Tail  Up 
Wing,  + Tip  Up 
Horizontal  Tail,  + Tip  Up 

Fuselage,  + Tail  Right 
Fin,  + Tip  Right 

Wing,  + Tip  Aft 

+ Nose  Up 

+ Left  Wing  Up 

+ Nose  Right 

+ Up 

+ Right* 

+ Leading  Edge  Up 
+ Trailing  Edge  Down  or  Left 
+ Nose  Up 
+ Nose  Left 


*Lateral  acceleration  at  BS  1655  is  + Left 
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SECTION  IV 


FLIGHT  TEST  CONDITIONS 

The  data  acquired  for  application  of  parameter  estimation  methods  was  obtained 
from  the  CCV  . light  Test  Program.  The  nominal  test  conditions  were;  an 
altitude  of  21,000  feet,  an  airspeed  of  305  KCAS,  and  a gross  weight  of 
260,731  pounds.  The  test  conditions  for  the  eight  flight  tests  selected  are 
shown  in  Table  2.  For  seven  of  the  eight  test  conditions,  the  Lateral  Aug- 
mental  Stability  (LAS)  system  was  in  operation.  For  the  rudder  transient 
input  condition  (45.2.5)  the  LAS  was  not  operating. 


TABLE  2.  FLIGHT  TEST  CONDITIONS 


t> 


SECTION  V 

FLIGHT  TEST  DATA  AND  ANALYTICAL  COMPARISONS 


This  section  presents  plots  o£  the  flight  test  data  selected  as  being  suitable 
for  parameter  estimation  methods;  with  plots  of  analytical  responses  calculated 
for  each  test  condition.  Flight  test  and  analytical  response  data  are  presented 
in  Figures  4 through  8 for  the  five  flight  test  conditions  with  transient  con- 
trol inputs.  Flight  test  data  for  the  three  flight  test  conditions  with 
frequency  sweep  control  Inputs  is  presented  in  Figures  9,  10  and  11.  Analytical 
frequency  response  calculations  were  made  for  each  of  the  frequency  aweep 
conditions . 

5.1  Flight  Test  Data 

The  plots  of  measured  flight  test  data  presented  in  this  section  are  for 
the  eight  flight  test  conditions  shown  in  Table  2.  The  time  history 
plots  of  the  measured  responses  represent  only  the  dynamic  part  of  the 
response  (i.e.,  the  initial  values  of  the  responses  were  subtracted  from 
the  respective  responses  prior  to  plotting).  The  plots  of  control  posi- 
tion represent  the  actual  control  position.  For  ease  of  viewing,  all 
of  the  plotted  test  data  were  smoothed  via  a three-point  running  average 
of  the  data. 

5.2  Analytical  Comparisons 

Analytical  responses  were  calculated  for  the  five  flight  conditions  with 
transient  control  inputs.  Responses  were  calculated  for  a period  of  ten 
seconds  for  each  condition;  even  though  the  flight  tests  varied  in  length 
from  5.2  to  16.4  seconds. 

The  analytical  responses  were  computed  from  the  equations  of  motion  and 
response  coefficients  presented  in  Appendix  B.  The  computations  were 
made  using  a linear  time  response  program  which  uses  a modified  Fast 
Fourier  Transform  solution.  For  the  symmetric  control  input  conditions, 
the  equations  of  motion  were  used  as  presented.  However,  the  anti- 
symmetric equations  of  motion  for  the  airplane  without  yaw  damper  did 
not  have  sufficient  Dutch  Roll  stability.  As  a simple  correction  for 

niirnnooc  nf  voyH  ferine*  t*V«A  fMoVit*  hoot*  r)ai*a  r*Ar1np  nn  nronooq  o 

— j o — — o— — 

tious  yaw-rate- to-rudder  damper  was  added.  The  yaw  damper  used  was  that 
used  on  the  CCV  airplane,  but  with  gain  reduced  by  two-thirds.  The 
damper  equation  is  shown,  in  a block  diagram  form,  in  Appendix  B. 

Responses  calculated  using  the  basic  symmetric  equations  of  motion  and 
the  augmented  antisymmetric  equations  of  motion  are  presented  with  the 
corresponaing  tp  ,t  data  in  Figures  4 through  8.  These  response  plots 
show  reasonable  agreement  between  the  flight  test  measurements  and  the 
analytical  responses. 

It  was  concluded  from  this  comparison  that  the  flight  test  data  reduc- 
tion process  was  satisfactory,  that  sensor  identification  was  correct, 
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and  that  the  calibration  ana  algebraic  sense  were  correct  for  all  re- 
sponses. Further,  it  was  concluded  that  the  airplane  structural  and 
geometric  description  given  in  the  appendices,  when  used  with  computer 
programs  such  as  FLEXSTAB  2 (which  predicts  aerodynamic  influence  and 
forms  equations  of  motion  and  response  equations)  is  satisfactory  for 
the  verification  of  such  programs. 

5.3  Analytical  Frequency  Responses 

Analytical  frequency  responses  were  calculated  using  the  control  inputs 
from  the  three  frequency  sweep  flight  test  conditions.  The  amplitudes 
used  for  the  control  surface  deflection  were  the  same  as  those  used  in  the 
flight  tests.  However,  for  the  frequencv  sweep  the  analytical  calculations 
were  started  at  .005  Hertz  rather  than  the  .5  Hertz  starting  frequency 
used  in  the  flight  tests.  The  change  in  starting  frequency  allows  the 
effects  of  the  rigid  body  modes  to  be  observed  in  the  frequency  response 
plots.  Frequency  responses  are  shown  in  Figures  12,  13  and  14  for  out- 
board aileron  input,  outboard  flaperon  input  and  vertical  canard  input, 
respectively. 

Each  of  the  frequency  response  plots  presents  two  items  of  data  as  a 
function  of  frequency.  The  data  referenced  to  the  left  hand  scale  is 
the  amplitude  of  the  response.  The  units  for  the  response  amplitude  are 
given  at  the  top  rf  each  plot  along  with  the  response  identification. 

The  phase  relationship  between  the  response  and  the  input  is  referenced 
to  tl  e right  hand  scale  with  the  phase  angle  reduced  to  fall  within  the 
limits  of  the  ±180  degrees. 

In  order  to  observe  the  cyclic  nature  of  the  phase  angle,  trace  the  path 
of  the  plotted  symbols  as  it  disappears  at  the  top  (bottom)  of  the  plot 
and  reappears  at  the  bottom  (top).  For  each  complete  cycle  across  the 
plot,  add  (subtract)  360  to  (from)  the  plotted  angle. 

As  a guide  for  reading  the  frequency  response  plots,  Table  3 on  page  145, 
has  been  included  to  show  the  amplitude  and  phase  angle  for  each  of  the 
plotted  responses  at  a frequency  of  1.0  Hertz. 
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Figure  4.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.7  (Continued) 
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Figiire  4.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.7  (Continued) 
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Figure  5.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.17  (Continued) 
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Figure  5.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.17  (Continued) 
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Figure  5.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.17  (Continued) 
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Figure  5.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.17  (Continued) 
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Figure  5.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.17  (Continued) 
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Figure  5.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.17  /r^-itinued) 
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Figure  5.  Flight  Test  and  Analytical  Response  Comp 
for  Flight  Condition  44.3.17  (Continued) 
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Figure  6.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.9  (Continued) 
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Figure  6.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.9  (Continued) 
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Figure  6.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.9  (Continued) 
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Figure  6.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.9  (Continued) 
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Figure  6.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.3.9  (Continued) 
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Figure  7.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.48.2  (Continued) 
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Figure  7.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.48.2  (Continued) 
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Figure  7.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.48.2  (Continued) 
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FLIGHT  TEST  DATA 


Figure  7.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.48.2  (Continued) 


<yniw»iii  )■>  r»Wn  irtnAami; 


r 


mpmmr** 


J wikffl 


60 

5 

U 

u 

A 

bo 


r—i 

CO 

a 

•H 

4J 

0) 


FLIGHT  CONDITION  ’44.2.48.2 


Figure  7.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.48.2  (Continued) 
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Figure  7.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.48.2  (Continued) 


Figure  7.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  44.2.48.2  (Concluded) 
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Figure  8.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  45.2.5  (Continued) 
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Figure  8.  Flight  Test  and  Analytical  Response  Comparisons 
for  Flight  Condition  45.2.5  (Continued) 
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44.1.1  (Continued) 


44.1.1  (Continued) 
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Figure  9.  Flight  Test  Response  Data  for  Flight  Condition 
44.1.1  (Continued) 
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Figure  9.  Flight  Test  Response  Data  for  Flight  Condition 
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Figure  10.  Flight  Test  Response  Data  for  Flight  Condition 
44.1.5  (Continued) 
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Figure  10.  Flight  Test  Response  Data  for  Flight  Condition 
44.1.5  (Continued) 
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Figure  10.  Flight  Test  Response  u<ica  for  Flight  Condition 
44.1.5  (Continued) 
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Figure  10.  Flight  Test  Response  Data  for  Flight  Condition 
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Figure  10.  Flight  Test  Response  Data  for  Flight  Condition 
44.1.5  (Continued) 
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Figure  11.  Flight  Test  Response  Data  for  Flight  Condition 
45.1.1  (Continued) 
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Figure  11.  Flight  Test  Response  Data  for  Flight  Condition 
45.1.1  (Continued) 
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Figure  11.  Flight  Test  Response  Data  for  Flight  Condition 
45.1.1  (Continued) 
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TABLE  3.  REFERENCE  VALUES  OF  AMPLITUDE  AND  PHASE  ANGLE  FROM  EACH 
FREQUENCY  RESPONSE  PLOT  - VALUES  ARE  FOR  A FREQUENCY  OF 
1.0  HERTZ  (Concluded) 
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Figure  12.  Frequency  Responses  - Outboard  Aileron  Frequency  Sweep, 

.005  to  5.0  Hz,  1.6  Degree  Amplitude  Sine  Wave  Input  (Continued) 
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Figure  13.  Frequency  Responses  - Outboard  Flaperon  Frequency  Sweep, 

.005  to  5.0  Hz,  2.0  Degree  Amplitude  Sine  Wave  Input  (Concluded) 
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PARAMETER  ESTIMATION  BY  CORRELATIONS  AND  LEAST  SQUARES 


i 


s 


} 


f 

L 


f 


L 

i 

f 


i 

t 


I ! 


f 


Parameter  estimation  (PE)  methods  are  categorized,  in  Reference  7 , as  (1)  Equa- 
tion Error  Methods,  (2)  Output  Error  Methods,  and  (3)  Advanced  Methods.  In 
general,  the  first  category  is  quick-look  methods,  unable  to  account  for  meas- 
urement or  input  (process)  noise.  The  second  can  account  for  measurement  noise, 
and  the  third  can  account  for  both  measurement  and  some  types  of  input  noise. 

The  PE  method  used  in  this  project  is  not  readily  categorized  in  this  manner. 

It  is  basically  a least-squares  method,  which  would  be  the  first  category. 
However,  the  .least  squares  process  is  not  applied  to  the  airplane  motions 
it  is  applied  to  the  cross-correlation  functions  of  the  airplane  state  variables. 
In  this  way  it  is  possible  to  eliminate  the  effects  of  measurement  noise  and 
input  noise  of  any  spectral  shape.  As  with,  any  cross-correlation  process,  the 
ability  to  resolve  narrow  band  spectral  peaks  is  increased  with  a long  maximum 
correlation  lag.  If  the  signal/noise  ratio  is  poor,  many  statistical  degrees 
of  freedom  (ratio  of  total  test  time  to  maximum  correlation  lag)  will  be 
needed.  If,  in  fact,  there  is  no  noise  or  unknown  input  in  the  data,  the  PE 
method  becomes  equivalent  to  a weighted  least-squares  method. 

A critical  element  of  this  estimation  if  the  recognition  that  the  aerodynamic 
displacement  and  rate  derivatives  will  vary  with  frequency.  To  permit  the 
observation  of  this  variation,  the  state  vector,  normally  [sq,  qj  was  augmented 
with  a quasi-integral  to  become  [sq,  q,  sq/(s+/S)]  . This  combination  of 
freedoms  has  been  found  to  adequately  represent  theoretical  unsteady  aerody- 
namics over  the  frequency  range  of  interest  (0-3  Kz)  for  the  B-52. 

The  computational  effort  is  increased  with  this  addition.  An  eight  degree-of- 
freedom  (flight  modes  plus  elastic  modes)  airplane  would,  for  example,  have 
24  state  variables.  In  order  to  reduce  computing  costs,  the  PE  method  has  been 
programmed  to  perform  the  correlations  and  estimation  for  an  equivalent  third 
order  set  of  difference  equations.  Conversion  to  differential  equations  is  then 
accomplished  via  eigensystem  analysis. 

6.1  Development  of  Parameter  Estimation  Method 

The  following  paragraphs  show  how  the  generalized  displacements  (elastic 
mode  displacements)  were  obtained  from  measured  responses,  the  develop- 
ment of  the  difference  equation  coefficients  and  the  conversion  to 
differential  equation  coefficients, 

6.1.1  Generalized  Displacements  from  Test  Data,  A Truncated  Mode  Approximation 

In  a theoretical  dynamic  analysis,  the  structural  motions  are  assumed  to 
be  a superposition  of  center-of-gravity  motion?  plus  a limited  number  of 
elastic  mode  deformations.  The  linear  and  angular  motions  of  points  on 
the  structure  can  be  written 

{d(t)f  = [P]  {q(t)f 
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where  |q(t)[  are  the  instantaneous  mode  amplitudes  and  the  columns  of 
[p]  are  deflections  per  unit  modal  displacement.  With  sensors  chosen  as 
accelerometers,  the  sensor  outputs  would  be  related  to  the  elastic  and 
center-of-gravity  motions  by 


jd(t)|  > [P]|q'(t)( 

where  |q(t)}  are  the  instantaneous  modal  accelerations. 


(2) 


If  strain  gages  provide  the  measured  response,  the  loads  are  related  to 
the  elastic  deformations  by 

|£(t)}  = [l]  {q(t)}  (3) 

where  the  columns  of  Jl]  are  mode  displacement  load  coefficients.  The 
columns  of  [p]  and  [lJ  are  determined  from  a theoretical  vibration 
analysis.  When  using  either  [p]  or  [l]  it  is  assumed  that  the  matrix 
will  contain  more  rows  than  columns.  Using  strain  gage  data,  the 
modal  displacements  can  be  defined  by  the  generalized  inverse  as  follows: 


(4) 


For  symmetric  tests,  the  state  variables  could  be  angle-of-attack, 
pitch  rate,  and  modal  displacements  for  six  elastic  modes.  Then  the 
state  vector  is 

a (t)\ 

« » * (t)l 

|x(t) J-  - qj(Y) 

(qj.(t) 


(5) 


For  the  antisymmetric  test,  the  variables  could  be  side  slip  angle,  roll 
rate,  yaw  rate  and  modal  displacements  for  five  elastic  modes. 

6.1.2  Estimation  of  Difference  Equation  Coefficients 

Coefficients  are  sought  for  the  third  order  difference  equation, 

■ Fi]  N + Fdh-il  + F3])*i-2f  (6> 

+Pi]  fd+Pd  + LD3j  iW  + pii 

Where  the  state  vector  |x|  has  been  calculated  from  measured  data  and 
the  inputs  If}  have  been  measured.  All  measurements  are  assumed  to  be 
noisy.  The  vector  |c|  contains  the  unmeasured  inputs  (process  noise). 

When  the  above  equation  is  rewritten  in  the  following  form, 


|xi+l(  = [C1  j C2  | C3  j D1  j °2  j °3]  / xi 


xi-l 

xi-2 

fi 

fi-1 

fi-2 


i€ii 
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the  following  identity  is  readily  established: 

I n ^ ® tv  1 n I 


[Uk+3]  - {pl  | C2  ; c3  | Dx  i D2  : D3] 


u, 


k+2 


u, 
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+ M 


(8) 


Where  IU],  [V],  and  [E]  are  cross-correlation  matrices  defined  as 
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The  measurement  noise  vanishes  in  the  correlation  calculation  if  n-»eo. 
The  matrix  [e]  will  vanish  relative  to  [u]  and  [v]  if  the  jc}  are 
uncorrelated  with  the  |f|  and  n/m->oo.  m must  be  at  least  large 
enough  to  provide  the  aesired  s-plane  resolution,  as  defined  by  the 
Nyquist  criteria.  The  number  of  statistical  degrees  of  freedom,  n/m, 
required  will  depend  upon  the  amount  of  noise  and  undesired  inputs 
present. 


The  difference  equation  coefficients  are  then  found  by  least  squares, 
namely : 

/m-3  „ \ /m-3  \ -1 
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(12) 


The  preferred  test  for  this  method  of  analysis  would  be  independent 
random  inputs  to  each  of  the  control  surfaces,  spectrally  she ped  to 
provide  uniform  signal/noise  over  the  test  bandpass. 

6-1.3  Conversion  to  Differential  Equation  Coefficients 

After  the  coefficient  matrices  [c]  and  [d]  of  the  difference  equations 
have  been  determined,  the  coefficient  matrices  [a]  and  [b]  of  the 
equivalent  differential  equations  may  be  calculated. 
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The  required  derivations  follow  from  a comparison  of  the  following  first 
order  forms: 
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where  q- 


±11 


— »r  and  f- 


+ 0 


Taking  the  LaPlace  transform  of  each  shows  the  equivalence  of  the  two 
eigensys terns: 
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Thus  AA  = ^-  L0Ge  (XC} 

.1  _ i 


and  [A1!A2JA3]  = [l;0,.0][T]  [diag  (XA)]  [t]”1,  Where  [t]  = 


and  where  jx^Jis  a partition  of  the  corresponding  vector 
of  the  C-matrix. 
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Although  [Aj,  [a£]  , and  [A3]  are  completely  determined  by  [C]]  , [03], 
[C3]  and  the  specification  of)S,  the  determination  of  [Bj , [B2],  and 
[B3]  requires  an  assumption  about  the  form  of  tne  forcing  function. 
For  example,  if 

W2 


4+1 
h-i] 

and,  using  g 


r 

2 “ 


0 

At 


”2”"  J 


4tl 


1 j. — 1 


-dtl 


s f 


s +fi 


<i  -f) 


0 

M 

CM 

j 

r 

! « 
1 

c 

0 

[° 

— « 
I 

0 

_ 

we  can  write 
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so  that 
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6.2  Parameter  Estimation  Data 


The  parameter  estimation  method  described  in  the  previous  paragraphs  was 
applied  to  the  three-frequency  sweep  and  the  five  transient  flight  test 
conditions  described  in  Section  IV.  The  eigenvalues  obtained  from  the 
frequency  sweep  test  data  compare  much  more  favorably  with  the  theoretical 
eigenvalues  than  the  ones  obtained  from  the  transient  test  data.  The 
differences  are  due  in  part  to  the  amount  of  data  available.  The  fre- 
quency sweep  test  were  135  to  140  seconds  in  duration,  while  the  transient 
test  lasted  from  5 to  16  seconds.  A discussion  of  the  results  of  the 
parameter  estimation  is  given  in  the  following  paragraphs. 

6.2.1  Estimated  Data  from  Outboard  Aileron  Sweep 

Application  of  parameter  estimation  to  the  test  involving  the  outboard 
ailerons  produces  the  most  satisfactory  results.  This  should  not  be 
entirely  unexpected  since  the  outboard  ailerons  excite  the  lower  fre- 
quency structural  modes  to  a greater  extent  than  either  the  outboard 
flaperons  or  the  vertical  canard.  A plot  of  theoretical  and  estimated 
system  eigenvalues  is  shown  in  Figure  15.  Since  the  control  input  fre- 
quency sweep  started  at  about  3 radians/sec,  the  short  period  mode  (1.62 
radians/sec)  was  not  well  defined.  The  candidate  modes  included  in  the 
parameter  estimation  included  the  first  six  elastic  modes  (6  to  16  radians/ 
sec).  Therefore,  eigenvalues  above  16  radians/sec  were  not  well  defined. 
The  estimated  eigenvalues  shown  in  Figure  15  were  derived  from  bending 
moment  measurements,  angle-of-attack,  and  center-of-gravity  pitch  rate. 
There  were  no  bending  moment  measurements  for  the  nacelle  struts  and,  as 
expected,  the  nacelle  modes  at  about  12  radians/sec  were  not  well  defined. 

Flight  Condition  44.2.48.2  was  transient  input  to  the  outboard  aileron. 

The  input  to  the  transient  condition  was  applied  to  the  system  coefficients 

derived  from  the  frequency  sweep  data,  resulting  in  an  approximation  of 

the  measured  data  for  the  transient  condition.  Figure  16  shows  the 
generalized  freedoms  calculated  using  the  derived  coefficients  along 
with  the  input  used.  The  computed  elastic  mode  displacements  were  used  to 
generate  the  bending  moment  time  histories  shown  in  Figure  17.  Con- 
sidering the  differences  between  the  theoretical  eigenvalues  and  the 

estimated  eigenvalues,  the  reconstructed  bending  moments  do  not  appear 
to  be  unreasonable. 

An  attempt  was  made  to  resolve  the  nacelle  modes  using  the  acceleration 
data,  since  both  vertical  and  lateral  accelerations  were  measured  for 
the  Inboard  and  outboard  nacelles.  The  estimated  system  eigenvalues 
shown  in  Figure  18  do  not  compare  very  favorably  with  the  theoretical 
eigenvalues  nor  with  the  eigenvalues  estimated  from  the  bending  moment 
measurements.  This  may  indicate  that  the  acceleration  measurements  are 
not  as  useful  in  resolving  elastic  mode  deformations  as  are  bending 
moments. 

6.2.2  Estimated  Data  from  Outboard  Flaperon  Frequency  Sweep 

The  eigenvalues  estimated  from  the  outboard  flaperon  frequency  sweep 
c I ! test  data  are  shown  in  Figure  19.  The  eigenvalues  were  estimated  using 
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bending  moment,  angle  of  attack  and  pitch  rate , The  estimated  eigen- 
values show  rather  poor  comparison  with  the  theoretical  eigenvalues. 
However,  the  flaperons  do  not  excite  the  lower  frequency  modes  very 
well  and  the  comments  about  the  short  period  and  higher  frequency  modes 
and  lack  of  nacelle  strut  data  still  apply.  Since  there  is  no  transient 
flaperon  input  tests  included  in  the  data  package,  no  responses  were 
computed  using  the  derived  system  coefficients. 

6.2.3  Estimated  Data  from  Vertical  Canard  Frequency  Sweep 

The  eigenvalues  estimated  from  the  vertical  canard  frequency  sweep  test 
data  are  shown  in  Figure  20.  For  this  antisymmetric  condition  the  data 
used  for  the  estimation  was  sideslip  angle,  roll  and  yaw  rates,  and  bend- 
ing moment  measurements.  The  bending  moment  data  was  used  to  compute 
elastic  mode  deflections  for  the  first  five  elastic  modes  (11  to  14.1 
radians/sec) . The  elastic  mode  deflections  were  then  used  with  the 
sideslip  and  roll  and  yaw  rates  for  the  parameter  estimation  routine. 

As  with  the  flaperon  data,  the  estimated  eigenvalues  do  not  compare 
very  well  with  the  theoretical  eigenvalues.  This  is  in  part  due  to  the 
lack  of  nacelle  data,  lack  of  input  below  the  roll  convergence  mode 
frequency  and  the  fact  that  the  vertical  canard  does  not  excite  the 
lower  frequency  elastic  modes  sufficiently  to  provide  good  resolution 
of  the  modes. 

6.2.4  Estimated  Data  from  the  Transient  Input  Test  Conditions 

Eigenvalues  were  estimated  from  the  recorded  data  from  each  of  the  tran- 
sient input  flight  tests  (elevator  sine  wave,  elevator  rectangle,  ele- 
vator triangle,  outboard  aileron  and  rudder  sine  wave) . The  estimated 
eigenvalues  are  shown,  along  wirh  the  theoretical  eigenvalues,  in  Figures 
21  through  25.  The  plots  of  estimated  and  theoretical  eigenvalues  show 
very  poor  comparison  between  the  two.  There  are  two  reasons  for  the  poor 
comparison.  One  reason  is  that  there  is  insufficient  data  available  from 
the  short  duration  transient  tests  to  eliminate  the  effects  of  noise 
(either  processing  or  external) . The  second  reason  for  the  poor  comparison 
is  that  the  transient  control  inputs  do  not  excite  all  of  the  rigid  body 
and  elastic  modes  being  considered  in  the  parameter  estimation  model. 
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Figure  15.  Estimated  Eigenvalues,  From  Outboard  Aileron  Frequency 
Sweep  Test  - Bending  Moment  Data 
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Figure  17.  Bending  Moments  for  Outboard  Aileron  Transient,  Derived  from 
System  Coefficients  Estimated  from  the  Outboard  Aileron 
Frequency  Sweep  Condition  - Bending  Moment  Data 
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Figure  17.  Bending  Moments  for  Outboard  Aileron  Transient,  Derived  from 
System  Coefficients  Estimated  from  the  Outboard  Aileron 
Frequency  Sweep  Condition  - Bending  Moment  Data  (Concluded) 
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Figure  18.  Estimated  Eigenvalues,  from  Outboard  Aileron  Frequency 
Sweep  Test  - Acceleration  Data 
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Figure  19.  Estimated  Eigenvalues,  from  Outboard  Flaperon  Frequency 
Sweep  Test  - Bending  Moment  Data 
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Figure  20.  Estimated  Eigenvalues,  from  Vertical  Canard  Frequency 
Sweep  Test  - Bending  Moment  Data 
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i Figure  21.  Estimated  Eigenvalues  from  Elevator  Sine  Wave  Input 

I Transient  Test  - Bending  Moment  Data 
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CCV  FLIGHT  TEST  DATA  - CONDITION  44.3.9  - ELEVATOR  TRIANGLE  TRANSIENT 
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Estimated  Eigenvalues  from  Elevator  Triangle  Input 
Transient  Test  - Bending  Moment  Data 
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Figure  24.  Estimated  Eigenvalues  from  Outboard  Aileron  Sine  Wave 
Input  Transient  Test  - Bending  Moment  Data 
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Figure  25.  Estimated  Eigenvalues  from  Rudder  Sine  Wave  Input 
Transient  Test  - Bending  Moment  Data 
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SECTION  VII 


CONCLUSIONS 


The  flight  test  and  analytical  data  documented  in  this  report  ar%  adequate  for 
parameter  estimation  applications  and  transient  response  analyses.  It  was 
concluded  from  the  comparisons  in  Section  V that  the  flight  test  data  reduction 
process  was  satisfactory.  The  order  of  magnitude  and  algebraic  sign  of  each 
response  agrees  with  analytical  predictions.  It  was  also  concluded  that  the 
structural  and  geometric  description  given  in  the  appendices  is  satisfactory 
for  the  verification  of  programs  that  predict  aerodynamic  influence  and  form 
equations  of  motion  and  response  equations. 

It  can  be  seen  from  the  plots  of  estimated  system  eigenvalues  that  the  flight 
test  data  does  indeed  contain  information  relative  to  several  structural  elastic 
modes.  However,  it  is  concluded  that  the  test  inputs  were  not  particularly 
well  suited  for  parameter  estimation,  since  the  individual  control  input  did 
not,  in  general,  excite  all  of  the  rigid  body  and  structural  modes  being 
considered . 
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SECTION  VIII 


RECOMMENDATIONS 


The  techniques  used  for  selecting  and  measuring  state  vectors,  for  permitting 
frequency  variations  in  aerodynamic  derivatives,  for  parameter  estimation  and 
the  use  of  finite  difference  equations  followed  by  conversion  to  differential 
equation  coefficients  are  recommended  to  future  investigators  of  parameter 
estimation  applied  to  flexible  aircraft. 

The  individual  control  input  conditions  from  the  CCV  Flight  Tests  did  not  tend 
to  excite  all  of  the  rigid  body  and  elastic  modes  under  consideration.  It  is 
recommended  that  tests  being  conducted  with  consideration  of  parameter  estima- 
tion application  be  conducted,  at  least  in  part,  with  random  uncorrelated 
multiple  inputs.  This  procedure  also  is  helpful  when  the  measured  data  are 
noisy  or  there  are  external  random  inputs  (i.e.,  atmospheric  turbulence).  The 
input  could  be  spectrally  shaped  to  yield  uniform  input  to  expected  noise 
ratios  over  the  bandpass  of  interest. 

For  parameter  estimation  applications  to  flexible  aircraft  test  data,  it  i_ 
recommended  that  structural  deformation  instrumentation  (i.e.,  strain  gages) 
be  a significant  part  of  the  instrumentation  package,  and  that  deformation 
data  be  collected  from  all  flexible  components  of  the  aircraft,  including 
nacelles. 

The  parameter  estimation  data  presented  in  this  document  is  not  sufficient  to 
define  a mathematical  model  but  is  used  to  illustrate  that  the  parameter  esti- 
mation method  presented  in  Section  VI  is  a viable  method.  For  the  method  to 
be  of  general  use  it  is  necessary  to  include  the  conversion  from  finite  differ- 
ence equation  coefficients  to  differential  equation  coefficients,  as  shown  in 
Paragraph  6.1.3. 
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APPENDIX  A 


FLEXIBLE  AIRPLANE  ANALYSIS  MODEL 


Introduction 

Appendix  A contains  a numerical  description  of  the  B-52  test  airplane 
structural  and  aerodynamic  configuration.  These  data  were  the  basis  of 
the  CCV  program  theoretical  analysis  and  system  design.  The  data  are 
presented  in  a form  which  will  be  readily  adaptable  to  most  flutter  and 
aero elastic  analysis  programs. 

Total  airplane  inertia  data,  results  of  the  structural  vibration  analy- 
sis and  inflight  eigenvalues  are  presented  in  Section  2.0.  These  data 
will  provide  a convenient  check  for  those  who  independently  model  the 
airplane.  Both  symmetric  and  antisymmetric  results  are  shown  for  the 
analysis  gross  weight. 

Section  3.0  contains  the  beam  element  description  which  represents  the 
elastic  structure.  Nodal  mass  data  are  given  for  the  empty  airplane 
and  separately  for  the  fuel  loading. 

Section  4.0  contains  the  aerodynamic  patch  description  used  for  the 
doublet  lattice  aerodynamic  analysis.  Since  the  origin  of  the  coordi- 
nate system  used  for  the  aerodynamic  patches  differs  from  that  of  the 
structural  analysis,  the  aerodynamic  coordinates  of  the  structural 
reference  points  used  to  transfer  motions  and  forces  are  also  given. 
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2.0  Total  Mass,  Component  and  Total  Vibration  Data 





TABLE  4.  TOTAL  AIRPLANE  MASS  AND  VIBRATION  DATA 


WEIGHT  = 260,731  LBS. 
BS  C.G.  = 863.74  In 
BWL  C.G.  = 189.86  In 

I = 4. 20  x 1010  LB-In2 
xx  iq  2 

i I =2.87  xl0lu  LB-In Z 

yy  in  o 

I = 6.99  x 101U  LB-Li 

ZZ  ft  9 

I = -7.46  X 10  LB-In 

Vff 


MODE 

FREQUENCY  ~ Hz 
SYM. 

A/S 

1 

.88 

1.36 

2 

1.74 

1.64 

3 

2.01 

1.89 

4 

2.15 

2.06 

5 

2.41 

2.29 

6 

2.50 

2.73 

7 

2.97 

2.92 

8 

3.14 

3.30 

9 

3.77 

3.69 

10 

5.78 

4.89 

11 

5.99 

5.38 

12 

6.07 

5.88 

13 

7.06 

5.92 

14 

7.59 

6.01 

15 

7.88 

6.52 

16 

10, 07 

7.66 

17 

10.81 

9.30 

18 

12.63 

10.63 

19 

13.48 

11.74 

20 

14.22 

12.05 

21 

16.05 

12.88 

22 

16.74 

13.95 

23 

17.03 

15.95 

24 

20,07 

16.16 

25 

21.34 

16.80 

26 

2L95 

17.87 

27 

24.18 

18.20 
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TABLE  5-  COMPONENT  VIBRATION  DATA 


FORE  BODY 

FREQUENCY-  Hz 

MODE 

SYM. 

A/S 

1 

4.63 

3.81 

2 

21.02 

12.39 

3 

AFT  BODY* 

46.66 

FREQUENCY  - Hz 

MODE 

SYM. 

A/S 

l 


i 


1 

1.35 

1.22 

2 

5.94 

2.92 

3 

11.39 

5.02 

4 

15.25 

8.38 

5 

27. 56 

14.46 

6 

37.28 

16.37 

7 

43.67 

21.20 

8 

VERTICAL  TAIL 

59.81 

FREQUENCY  - Hz 

MODE 

SYM. 

A/S 

1 

- 

5.16 

2 

- 

14.41 

3 

- 

17.93 

4 

WING** 

FREQUENCY  - Hz 

27.39 

MODE 

SYM. 

A/S 

1 

.75 

.78 

2 

V 

1.78 

1.79 

3 

2.20 

2.21 

4 

2.50 

2.55 

5 

3.57 

3.57 

n 

o 

4.30 

4.42 

7 

6.53 

6.57 

8 

6.85 

6.85 

9 

9.84 

9.95 

10 

12.36 

- 

11 

14.59 

- 

12 

16.80 

- 

13 

30.49 

- 

14 

21.23 

- 

15 

24.67 

- 

*Aft  body  component  modes  include  mass  of  rigid  H.  Tail  and  V.  Tail. 
**Wing  component  modes  include  mass  of  rigid  inboard  and  outboard  nacelle 
and  tank. 
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i TABLE  5.  COMPONENT  VIBRATION  DATA  (Concluded) 


HORIZONTAL  TAIL 

FREQUENCY  ~ Hz 

i MODE 

SYM. 

A/S 

1 

7.22 

7.22 

j 

2 

16.44 

16.44 

3 

21.01 

21.01 

j 

1 

4 

INBOARD  NACELLE 

32.27 

FREQUENCY  - Hz 

MODE 

SYM. 

A/S_ 

1 

2,04 

2.04 

2 

4.07 

4.07 

j 

3 

OUTBOARD  NACELLE 

6.01 

FREQUENCY  ~ Hz 

6.01 

1 

j MODE 

SYM. 

A/S 

1 

1 

2.10 

2.10 

J 

2 

4.02 

4.02 

3 

1 

1 

i 

5.94 

5.94 

i 

! TABLE  6.  THEORETICAL  SYSTEM  EIGENVALUES  - LIGHT  GROSS  WEIGHT, 

i i 21,000  FT.,  305  KCAS,  SYMMETRIC 


Real 

Imaginary 

Frequency 

Percent  of 

Root 

Part 

Part 

Hertz 

Critical  Damping 

1 

Short  Period  Mode 

-1.272 

1.618 

.258 

61.81 

Elastic  Mode  1 

-0.633 

6.090 

.969 

10.33 

' 

Elastic  Mode  2 

-0.384 

12.190 

1.940 

3.15 

Elastic  Mode  3 

-0.133 

12.556 

1.998 

1.06 

i 

Elastic  Mode  4 

-1.140 

14.531 

2.313 

7.82 

! 

Elastic  Mode  5 

-0.136 

14.634 

2.329 

0.93 

Elastic  Mode  6 

-0.340 

15.723 

2.502 

2.16 

i 

Elastic  Mode  7 

-0.926 

19.168 

3.051 

4.32 

1 

Elastic  Mode  8 

-0.216 

19.600 

3.119 

1.10 

Elastic  Mode  9 

-1.146 

23.954 

3.812 

4.78 

| 

Elastic  Mode  10 

-0.510 

36.387 

5.791 

1.40 

Elastic  Mode  11 

-0.213 

37.630 

5.989 

0.57 

1 

Elastic  Mode  12 

-0.949 

38.008 

6.049 

2.49 

1 

Elastic  Mode  13 

-1.018 

44.379 

7.063 

2.29 

! Elastic  Mode  14 

-5.793 

49.340 

7.853 

11.66 

Elastic  Mode  15 

-0.703 

49.346 

7.854 

1.42 

Elastic  Mode  16 

-1.862 

62.803 

9.996 

2.96 

Elastic  Mode  17 

-1.342 

67.468 

10.738 

1.99 

Elastic  Mode  18 

-1.209 

78.619 

12.513 

1.54 

Elastic  Mode  19 

-0.861 

84.463 

13.443 

1.02 

Elastic  Mode  20 

-2.300 

89.082 

14.178 

2.58 

Elastic  Mode  21 

-1.883 

98.862 

15.734 

1.90 

Elastic  Mode  22 

-1.672 

104.321 

16.603 

1.60 

Elastic  Mode  23 

-1.519 

106.223 

16.906 

1.43 

Elastic  Mode  24 

-1.974 

125.700 

20.006 

1.57 

Elastic  Mode  25 

-0.991 

130.777 

20.814 

0.76 

Elastic  Mode  26 

-2.025 

135.579 

21.578 

1.49 

Elastic  Mode  27 

-1.476 

151.173 

24.060 

0.98 

1 

| Roots  resulting  from  lift-growth 

terms  not 

shown. 
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TABLE  7.  THEORETICAL  SYSTEM  EIGENVALUES  - LIGHT  GROSS  WEIGHT 
21,000  FT.,  305  KCAS,  ANTISYMMETRIC 


Real 

Imaginary 

Frequency 

Percent  of 

Root 

Part 

Part 

Hertz 

Critical  Damping 

Roll  Convergence  Mode 

-1.100 

0.0 

0.0 

100.0 

Dutch  Roll  Mode 

-0.113 

1.342 

0.214 

8.36 

Elastic  Mode  1 

-0.551 

9.542 

1.519 

5.77 

Elastic  Mode  2 

-0.464 

11.064 

1.761 

4.19 

Elastic  Mode  3 

-0.090 

12.794 

2.036 

0.70 

Elastic  Mode  4 

-1.454 

13.548 

2.156 

10.67 

Elastic  Mode  5 

-0.247 

14.074 

2.240 

1.75 

Elastic  Mode  6 

-0.182 

17.887 

2.847 

1.02 

Elastic  Mode  7 

-3.276 

18.052 

2.873 

17.85 

Elastic  Mode  8 

-0.998 

21.192 

3.372 

4.71 

Elastic  Mode  9 

-0.759 

23.012 

3.662 

3.30 

Elastic  Mode  10 

-0.925 

30.502 

4.855 

3.03 

Elastic  Mode  11 

-1.095 

33.741 

5.370 

3.24 

Elastic  Mode  12 

-0.230 

37.003 

5.889 

0.62 

Elastic  Mode  13 

-0.252 

37.473 

5.964 

0.67 

Elastic  Mode  14 

-5.736 

40.637 

6.468 

13.98 

Elastic  Mode  15 

-0.389 

40.919 

6.513 

0.95 

Elastic  Mode  16 

-2.195 

48.096 

7.655 

4.56 

Elastic  Mode  17 

-1.263 

58.145 

9.254 

2.17 

Elastic  Mode  18 

-2.049 

66.201 

10.536 

3.09 

Elastic  Mode  19 

-1.338 

73.351 

11.674 

1.82 

Elastic  Mode  20 

-1.707 

75.182 

11.966 

2.27 

Elastic  Mode  ?l 

-1.259 

80.289 

12.778 

1.57 

Elastic  Mode  22 

-2.278 

87.165 

13.873 

2.61 

Elastic  Mode  23 

-2.008 

98.325 

15.649 

2.04 

Elastic  Mode  24 

-2.440 

100.424 

15.983 

2.43 

Elastic  Mode  25 

-1.681 

102.928 

16.382 

1.63 

Elastic  Mode  26 

-1.013 

106.940 

17.020 

0.95 

Elastic  Mode  27 

-1.212 

111.648 

17.769 

1.09 

Roots  resulting  from  lift-growth  terms  are  not  shown. 
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3.0  Structural  Segment  Model,  Mass  and  Stiffness  Data 


The  CCV  airplane  structure  was  modeled  with  bar  elements  connected  at 
nodes.  The  weights  were  panelled  so  that  each  panel's  weight  could  be 
attached  to  the  structure  at  a node.  This  section  contains  a description 
of  the  nodes,  the  stiffness  of  the  bar  elements,  and  the  nodaL  mass 
properties. 


I 
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TABLE  8.  DEFINITION  OF  MASS  DATA  ORIENTATION  SYSTEMS 


TABLE  10.  AFT  BODY  WEIGHT  AND  STIFFNESS  DATA  (NO  USABLE  FUEL) 


1655  203.06  2535  -14.3  -32. S 4.12  3.54  1.68 


1791.96  286. 13  203.5 


1820.52  I 615. 10  38  .35  , U .000  .02-1  .020  . 008 


560  14.5  -10.50 


4.0  Aerodynamic  Patch  Model  and  Data 

Aerodynamic  lifting  surface  theory  was  incorporated  using  a doublet 
lattice  finite  element  model  to  stabilize  the  collocation.  The  aero- 
dynamic elements  (patches)  are  shown  on  pages  233  through  236. 

Each  patch  was  assigned  to  a structural  reference  point  from  which  it 
obtained  its  motion  and  to  which  it  transferred  itc  airloads.  Each 
patch  was  also  assigned  to  a gust  reference  point  from  which  it  obtained 
its  gust  boundary  condition  (with  a time  delay) . The  location  and 
orientation  of  the  structural  reference  points  and  the  location  of  the 
gust  reference  points  are  shown  on  pages  237  and  238. 
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TABLE  20.  AERODYNAMIC  PATCH  NUMERICAL  DESCRIPTION 
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TABLE  20.  AERODYNAMIC  PATCH  NUMERICAL  DESCRIPTION  (Continued) 
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TABLE  21.  AUXILIARY  DATA  FOR  AERODYNAMIC  PATCH  DESCRIPTION 


TURAL  REF.  STA. 

X 

Y 

Z 

1 

172.000 

0.0 

80.000 

2 

331.000 

0.0 

80.000 

3 

510.000 

0.0 

80.000 

4 

695.000 

0.0 

0.0 

5 

805.000 

0.0 

0.0 

6 

917.000 

0.0 

0.0 

7 

1135.000 

0.0 

0.0 

8 

1376.000 

0.0 

0.0 

9 

1581.000 

0.0 

0.0 

10 

1719.000 

0.0 

0.0 

11 

1684.  800 

94.910 

0.0 

12 

1733.300 

181.450 

0.0 

13 

1763.600 

235.540 

0.0 

14 

632.000 

410.000 

48.000 

15 

864.000 

720.000 

41.000 

16 

645.800 

102.700 

0.0 

17 

709.400 

205.400 

0.0 

18 

773.100 

308.000 

0.0 

19 

819.700 

383.100 

0.0 

20 

838.100 

410.000 

0.0 

21 

891.000 

487.500 

0.0 

22 

944.000 

565.000 

0.0 

23 

997.000 

642.500 

0.0 

24 

1049.900 

720.000 

0.0 

25 

1120.000 

822.500 

0.0 

26 

1190.000 

925.000 

0.0 

27 

1253.300 

1017.500 

0.0 

28 

1804.800 

94.910 

0.0 

29 

1825.000 

205.000 

0.0 

30 

uiu.uuu 

230. 000 

0.0 

31 

980.000 

430.000 

0.0 

32 

1185.000 

780.000 

0.0 

33 

1080.000 

700.000 

0.0 

34 

1656.800 

0.0 

-58. 100 

35 

1692.350 

0.0 

-134.900 

36 

1736.960 

0.0 

-231.200 

37 

1796.030 

0.0 

-358.730 

38 

1780. 100 

0.0 

-114.670 

39 

1810.700 

0.0 

-227.800 

40 

1841.300 

0.0 

-340.930 

ORIENTATION 
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STRUCTURAL  REFERENCE  POINT  ORIENTATION  SYSTEMS 
2,  3)  = (X,  Y,  Z)-  ANGLES  IN  DEGREES 


APPENDIX  B 


TEST  VEHICLE  THEORETICAL  EQUATIONS  OF  MOTION 


1.0  Introduction 


Appendix  B contains  a description  of  the  CCV  airplane  mathematical  model 
and  listings  of  the  coefficient  matrices  for  the  light-gross-weight  con- 
figuration. Included  in  the  mathematical  model  description  is  a block 
diagram  for  the  Electronic  Yaw  Damper  (EYD)  which  is  required  for  system 
stability  in  the  antisymmetric  analyses. 


2.0  Equations  of  Motion 

The  structural  equations  of  motion  (EOM)  were  furnished  in  the  following 
standardized  form. 


|s^  [Mass]  +s  [Damping]  + [Stiffness]  j |q(s)} 

+ (S2  P[CJ  + spu0  [c2]  +pu02  [c3]  +P0o2  £ [Dj 

+ h poi + "uo  £ m -»£*)  • ° 


Nondimensionalizing  the  equations  in  terms  of  p and  UQ  allows  a given 
tape  to  be  used  over  a moderate  range  of  airspeeds  and  altitudes.  The 
terms  in  s/(s+U0d^)  represent  unsteady  aerodynamic  effects.  Similar 
effects  are  modeled  for  gust  Influences.  Ordinarily,  there  are  30 
independent  coordinates  in  the  |q(s)|  plus  surface  positions.  Gust 
inputs  are  usually  modeled  for  nine  lift  panels.  The  independent 
coordinates  are  essentially  arbitrary  in  definition  except  for  the 
first  three,  which  are  the  rigid  body  coordinates,  lateral  or  longitudi- 
nal. The  others  are  dominant  structural  modes. 

Any  physical  motion  or  structural  load  is  a weighted  summation  of  the 
q's.  Therefore,  included  with  the  listings  of  EOM  coefficient  matrices 
is  a "modal  influence  matrix"  (symbol  (f) ) . Each  row  of  the  modal 
influence  matrix  represents  the  contributions  of  the  q's  to  the  motion 
of  or  the  load  at  a particular  point  on  the  airframe 

A block  diagram  of  the  Electronic  Yaw  Damper  (EYD)  is  shown  on  page  241. 
The  EYD  is  required  in  the  antisymmetric  analyses  in  order  to  have  a 
stable  dutch  roll  mode. 

The  coefficient  matrices  listed  on  the  following  pages  are  for  both  the 
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symmetric  and  antisymmetric  mathematical  models.  Each  of  the  mathemati- 
cal models  contain  three  rigid  body  degrees  of  freedom  and  27  elastic 
mode  degrees  of  freedom.  The  remaining  degrees  of  freedom  are  for  the 
control  surfaces.  In  the  symmetric  model  the  control  freedoms,  in  order, 
are  elevator,  outboard  aileron  and  outboard  flaperon.  For  the  anti- 
symmetric model  the  control  freedoms  are  rudder  and  vertical  canard. 
Summary  lists  of  the  matrices  printed  on  the  following  pages  are  shown 
on  pages  244  and  353.  The  arrangement  of  the  response  and  load  coeffi- 
cients in  the  matrix  (PHI)  is  shown  on  page  242. 

The  matrices  PUFF  through  FHON  are  arranged  so  that  each  column  repre- 
sents a generalized  freedom;  starting  with  the  rigid  body  freedoms 
(X,  Z and  -©•  for  the  symmetric  case  and  Y,  <p  and  ip  for  the  antisymmetric 
case)  followed  by  27  elastic  modes.  The  rows  of  each  matrix  represent 
modal  displacements  and  rotations.  The  rows  are  arranged  so  that  all  of 
the  X displacements  are  grouped  together  followed  by  all  of  the  Y dis- 
placements and  so  on  through  the  Z displacements  and  the(p,  -0-and  ip  rota- 
tions. For  the  three  degree-of-freedom  components  (forward  fuselage,  aft 
fuselage  and  vertical  tail)  the  null  displacements  and  rotations  have 
been  deleted  (Y,  <j>  and  \ p for  the  symmetric  case  and  X,  Z and  for  the 
antisymmetric  case) . 

Rigid  body  stability  derivatives  have  oeen  extracted  from  the  equation  of 
motion  matrices.  The  derivatives,  along  with  the  matrix  elements  from 
which  they  were  derived,  are  shown  on  pages  245  and  354  for  the  symmetric 
and  antisymmetric  models,  respectively. 


FIGURE  28.  ELECTRONIC  YAW  DAMPER  FUNCTIONAL  BLOCK  DIAGRAM 


TABLE  24 


RECORDED  FLIGHT  TEST  MEASUREMENTS 


1 


If 
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l. 
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i 


i 
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i 
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NO. 

ITEM 

LOCATION 

UNITS 

1 

Vertical  Bending  Moment 

BS760 

In. -LBS 

2 

Vertical  Bending  Moment 

BS1222 

) 

k 

3 

Vertical  Bending  Moment 

BS1412 

4 

Lateral  Bending  Moment 

BS760 

5 

Lateral  Bending  Moment 

BS1222 

6 

Lateral  Bending  Moment 

BS1412 

7 

Vertical  Bending  Moment 

LWS222 

8 

Vertical  Bending  Moment 

LWS820 

9 

Vertical  Bending  Moment 

LWS974 

10 

Chordwise  Bending  Moment 

LWS222 

11 

Chordwise  Bending  Moment 

LWS820 

12 

Chordwise  Bending  Moment 

LWS974 

13 

Vertical  Bending  Moment 

RWS222 

14 

Vertical  Bending  Moment 

RWS820 

15 

Vertical  Bending  Moment 

RWS974 

16 

Chordwise  Bending  Moment 

RWS222 

17 

Chordwise  Bending  Moment 

RWS820 

18 

Chordwise  Bending  Moment 

RWS974 

19 

Vertical  Bending  Moment 

LHTBL56 

20 

Vertical  Bending  Moment 

RHTBL56 

f 

21 

Lateral  Bending  Moment 

FS135 

I 

l.-LBS 

22 

Pitch  Rate 

BS860 

Deg/Sec 

23 

Roll  Rate 

BS860 

Deg/Sec 

24 

Yaw  R-*-e 

BS860 

Deg/Sec 

25 

Vertical  Acceleration 

BS'172 

g 

’s 

26 

Vertical  Acceleration 

BS860 

27 

Vertical  Acceleration 

BS1655 

28 

Lateral  Acceleration 

BS172 

29 

Lateral  Acceleration 

BS860 

30 

Lateral  Acceleration 

BS1655 

31 

Vertical  Acceleration 

LWBL565 

32 

Vertical  Acceleration 

LWBL925 

33 

Vertical  Acceleration 

LWS1359 

34 

Vertical  Acceleration 

Left  Ext.  Tank,  Nose 

35 

Vertical  Acceleration 

RW3L565 

36 

Vertical  Acceleration 

RWBL925 

37 

Vertical  Acceleration 

RWS1359 

38 

Vertical  Acceleration 

Right  Ext.  Tank,  Nose 

39 

Vertical  Acceleration 

Left  Inboard  Nacelle 

40 

Lateral  Acceleration 

Left  Inboard  Nacelle 

41 

Vertical  Acceleration 

Left  Outboard  Nacelle 

t 

42 

Lateral  Acceleration 

Left  Outboard  Nacelle 

g 

•a 

43 

Angular  Acceleration,  Y-Axis 

LWS540 

Rad/Sec^ 

44 

Angular  Acceleration,  Y-Axis 

LWS1359 

Rad/Sec^ 

45 

Angular  Acceleration,  Y-Axis 

RWS1359 

Rad/Sec^ 

46 

Delta  Inboard  Aileron 

- 

'Deg 

47 

Delta  Rudder 

- 

Deg 

i 
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TABLE  24.  RECORDED  FLIGHT  TEST  MEASUREMENTS  (Concluded) 


NO. 

ITEM 

LOCATION 

UNITS 

48 

Delta  Elevator 

— 

Deg 

49 

Delta  Vertical  Canard 

- 

J 

50 

Delta  Left  Horizontal  Canard 

— 

51 

Delta  Right  Horizontal  Canard 

- 

52 

Delta  Outboard  Aileron 

- 

53 

Delta  Outboard  Flaperon 

- 

54 

Delta  Inboard  Flaperon 

- 

55 

Angle  of  Attack 

Gust  Boom 

> 

' 

56 

Angle  of  Side  Slip 

Gust  Boom 

D< 

57 

Forward  Velocity 

Gust  Boom 

FT/Sec 

MEASUREMENT  SIGN  CONVENTION 


Vertical  Bending  Moments 


Lateral  Bending  Moment 

Chordwise  Bending  Moment 
Pitch  Rate 
Roll  Rate 
Yaw  Rate 

Vertical  Acceleration 
Lateral  Acceleration 
Angular  Acceleration,  Y-Axis 
Control  Surface  Deflection 
Angle  of  Attack 
Angle  of  Side  Slip 


Fuselage,  + Tail  Up 
Wing,  + Tip  Up 
Horizontal  Tail,  + Tip'Up 
Fuselage,  + Tail  Right 
Fin,  + Tip  Right 
Wing,  + Tip  Aft 
+ Nose  Up 
+ Left  Wing  Up 
+ Nose  Right 
+ Up 

+ Right* 

+ Leading  Edge  Up 
+ Trailing  Edge  Down  or  Left 
+ Nose  Up 
+ Nose  Left 


*Lateral  acceleration  at  BS  1655  is  + left 


I 

I 


I 


| 


i 
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TABLE  25.  EQUATION  OF  MOTION  AND  MODAL  COEFFICIENT  MATRICES  FOR 
SYMMETRIC  ANALYSES 


Matrix 

No.  of 

No.  of 

Name 

Rows 

Columns 

Description 

MASS 

33 

33 

Generalized  Mass  - Nominal  Gr.  Wt.,  260,731  Lbs. 

MASS 

33 

33 

Generalized  Mass  - Condition  44.2.7,  250,000  Lbs. 

MASS 

33 

33 

Generalized  Mass  - Condition  44.3.17,  258,900  Lbs. 

MASS 

33 

33 

Generalized  Mass  - Condition  44.3.9,  248,000  Lbs. 

MASS 

33 

33 

Generalized  Mass  - Condition  44.2.48.2,  242,000  Lbs 

MASS 

33 

33 

Generalized  Mass  - Condition  44.1.1,  272,300  Lbs. 

MASS 

33 

33 

Generalized  Mass  - Condition  44.1.5,  262,400  Lbs. 

Cl 

33 

33 

Aerodynamic  Influence  Coefficients  (acceleration) 

DAMP 

33 

33 

Generalized  Damping 

C2 

33 

33 

Aerodynamic  Influence  Coefficients  (velocity) 

STIF 

33 

33 

Generalized  Stiffness 

C3 

33 

33 

Aerodynamic  Influence  Coefficients  (displacement) 

D1 

33 

33 

First  Aerodynamic  Lag  Coefficients 

B1 

1 

33 

First  Aerodynamic  Lag  Constants 

D2 

33 

33 

Second  Aerodynamic  Lag  Coefficients 

B2 

1 

33 

Second  Aerodynamic  Lag  Constants 

RO 

33 

20* 

Gust  Influence  Coefficients 

R1 

33 

20* 

First  Gu&t  Lag  Coefficients 

BG1 

1 

20* 

First  Gust  Lag  Constant 

R2 

33 

20* 

Second  Gust  Lag  Coefficients 

BG2 

1 

20* 

Second  Gust  Lag  Constants 

R3 

33 

20* 

Third  Gust  Lag  Coefficients 

BG3 

1 

20* 

Third  Gust  Lag  Constants 

R4 

33 

20* 

Fourth  Gust  Lag  Coefficients 

BG4 

1 

20* 

Fourth  Gust  Lag  Constants 

XGST 

1 

10 

X Locations  of  the  Gust  Reference  Stations 

YGST 

1 

10 

Y Locations  of  the  Gust  Reference  Stations 

ZGST 

i 

10 

Z Locations  of  the  Gust  Reference  Stations 

PHI 

57 

30 

Modal  Coefficients  for  Flight  Test  Sensor  Stations 

PHFF 

24 

30 

Forward  Fuselage  Modal  Displacements  and  Slopes 

PHAF 

48 

30 

Aft  Fuselage  Modal  Displacements  and  Slopes 

PHWG 

96 

30 

Wing  Modal  Displacements  and  Slopes 

PHHT 

48 

30 

Horizontal  Tail  Modal  Displacements  and  Slopes 

PHIN 

12 

30 

Inboard  Nacelle  Modal  Displacements  and  Slopes 

PHON 

12 

30 

Outboard  Nacelle  Modal  Displacements  and  Slopes 

*This  matrix  is  partitioned  by  columns,  the  first  ten  columns  are  vertical 
gust  terms  and  the  last  ten  columns  are  lateral  gust  terms. 
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TABLE  26.  RIGID  AIRPLANE  STABILITY  DERIVATIVES  - SYMMETRIC  MODEL, 
MACH  = .687,  ALTITUDE  = 21,000  FT.,  CG  = .33  MAC 


Derivative  Equation  of  Motion  Coefficients 


Cr\ 

Da 

.001175/Deg 

2 

57. 3S 

Wl.3 

CL 

La 

.08736/'Deg 

- 2 
57. 3S 

W ,3 

c„ 

ma 

-.02198/Deg 

- 2 

Tn  1 

5.\3SC 

L 3J  3,3 

CLq 

6.4565/Rad 

- 4 f 
SC  [c2 

+ Tx  Dl  + 

1 

b2 

D*]  2,3 

Cmq 

-22.2672/Rad 

SC2  [C2 

+ Ti  D l + 

1 

B2 

d21 

3,3 

Ct 

L«e 

.003913/Deg 

57?3S  [ 

°3J  2,31 

c 

m«5e 

-.01194/Deg 

2 

57.3SC 

[°3]  3,31 

Ct 

<5  OA 

.001374/Deg 

572.SS  [ 

°3]  2,32 

C 

OA 

-.001361/Deg 

2 

57.3SC 

[°3]  3,32 

CT 

l<5  of 

.001602/Deg 

2 

57. 3S 

kl 

L •3-l  2,33 

Cm<5  OF 

.001246/Deg 

2 r< 

57.3SC  |_ 

=,] 

3,33 

S = 576 

2 

,000  IN.  = wing  area 

C = 275 

.5  IN.  = mean 

aerodynamic  chord 
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matrix  •mass*'  33  BY  33  SYMMETRIC  MATH  MODEL  NOMINAL  GROSS  WEIGHT,  260,731  LBS, 
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MATRIX  *MASS*  33  8V  33  SYMMETRIC  MATH  MODEL  FLIGHT  CONDITION  44.3.17.  258,900  UBS, 
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MATRIX  *MASS*  33  BY  33  SYMMETRIC  MATH  MODEL  FLIGHT  CONDITION  44.3.17,  258,900  LBS, 
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MATRIX  -MASS*  33  BY  33  SYMMETRIC  MATH  MODEL  FLIGHT  CONDITION  44.1.5,  262,400  LBS, 
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TABLE  27.  EQUATION  OF  MOTION  AND  MODAL  COEFFICIENT  MATRICES  FOR 
ANTISYMMETRIC  ANALYSES 


MATRIX 

NO.  OF 

NO.  OF 

NAME 

ROWS 

COLUMNS 

DESCRIPTION 

MASS 

32 

32 

Generalized  Mass  r Nominal  Gross  Weight  260,731  Lbs. 

MASS 

32 

32 

Generalized  Mass  - Condition  45.2.5  246,500  Lbs. 

MASS 

32 

'32 

Generalized  Mass  - Condition  45.1.1  254,000  Lbs. 

Cl 

32 

32 

Aerodynamic  Influence  Coefficients  (acceleration) 

DAMP 

32 

32 

Generalized  Damping 

C2 

32 

32 

Aerodynamic  Influence  Coefficients  (velocity) 

STIF 

32 

32 

Generalized  Stiffness 

C3 

32 

32 

Aerodynamic  Influence  Coefficients  (displacement) 

D1 

32 

32 

First  Aerodynamic  Lag  Coefficients 

B1 

1 

32 

First  Aerodynamic  Lag  Constants 

D2 

32 

32 

Second  Aerodynamic  Lag  Coefficients 

B2 

1 

32 

Second  Aerodynamic  Lag  Constants 

RO 

32 

20* 

Gust  Influence  Coefficients 

R1 

32 

20* 

First  Gust  Lag  Coefficients 

BG1 

I 

20* 

First  Gust  Lag  Constants 

R2 

32 

20* 

Second  Gust  Lag  Coefficients 

BG2 

1 

20* 

Second  Gust  Lag  Constants 

R3 

32 

20* 

Third  Gust  Lag  Coefficients 

BG3 

1 

20* 

Third  Gust  Lag  Constants 

R4 

32 

20* 

Fourth  Gust  Lag  Coefficients 

BG4 

1 

20* 

Fourth  Gust  Lag  Constants 

XGST 

1 

10 

X Locations  of  the  Gust  Reference  Stations 

YGST 

1 

10 

Y Locations  of  the  Gust  Reference  Stations 

ZGST 

1 

10 

Z Locations  of  the  Gust  Reference  Stations 

PHI 

57 

30 

Modal  Coefficients  for  Flight  Test  Sensor  Stations 

PHFF 

24 

30 

Forward  Fuselage  Modal  Displacements  and  Slopes 

PHAF 

48 

30 

Aft  Fuselage  Modal  Displacements  and  Slopes 

PHVT 

30 

30 

Vertical  Tail  Modal  Displacements  and  Slopes 

PHWG 

96 

30 

Wing  Modal  Displacements  and  Slopes 

PHHT 

48 

30 

Horizontal  Tail  Modal  Displacements  and  Slopes 

PHIN 

12 

30 

Inboard  Nacelle  Modal  Displacements  and  Slopes 

PHON 

12 

30 

Outboard  Nacelle  Modal  Displacements  and  Slopes 

*This 

matrix  is 

partitioned  by  columns,  the  first  ten  columns  are  vertical 

gust 

terms  and 

the  last  ten 

columns  are  lateral  gust  terms. 
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TABLE  28.  RIGID  AIRPLANE  STABILITY  DERIVATIVES  - ANTISYMMETRIC  MODEL, 


MACH  = .687, 

ALTITUDE  = 21,000  FT.,  CG  = . 

33  MAC 

Derivative 

Equation  of  Motion 

Coefficients 

cy 

-.09870/Deg 

57-3s  b]  i,3 

% 

-.001605/Deg 

57'3Sb  [°3  1 2,3 

c 

.002506/Deg 

57'3Sb  [°3]  3,3 

CyA 

yP 

-.01381/Rad 

ST  [ C2  + "B^  D1 

+ 4-  d2  1 

2 2J  1>2 

CpA 

*P 

-.4754/Rad 

- 4 r 1 

sb2  [ °2  + Bj  D1 

^ B ®2  1 

a2  1 J 2,2 

C A 
nP 

.02088/Rad 

’?"[ C 2 + V* 

+ V2  ]3,2 

c A 
Yr 

.4063/Rad 

sT  [C2  + D1 

*^]l.. 

.06667/Rad 

i7[°2  + V1 

+ k ”2  ] 2,3 

CnA 

-.1392/Rad 

^ h + k Di 

+ T„  D2  1 . , 

**  J Oy  J 

V 

.002207/Deg 

57'3S  [C3  ] 1,31 

cu 

wr 

.000185/Deg 

57-3Sb[  °3]  2,31 

c. 

r.  ; 

» r 

-.000881/Deg 

57-3Sb  [ C3  ] 3,31 
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TABLE  28.  RIGID  AIRPLANE  STABILITY  DERIVATIVES  - ANTISYMMETRIC  MODEL 

MACH  = .687,  ALTITUDE  = 21,000  FT.,  CG  = .33  MAC  (Concluded) 


Derivative 


Equation  of  Motion  Coefficients 


y<5 


vc 


vc 


n<5 


.000151/Deg 

-.000007/Deg 

.000047/Deg 


vc 


57. 3S 


57.3Sb 


57.3Sb 


b 

M,. 


32 


32 


S = 576,000  IN?  = wing  area 


b = 2220  IN.  = wingspan 
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APPENDIX  C 
MAGNETIC  TAPE  DATA 


1.0  Introduction 

Appendix  C presents  a description  of  the  data  tapes  containing  the 
B-52E  CCV  Flight  Test  Data  for  parameter  estimation  applications. 

2.0  Tape  Description 

The  magnetic  tapes  being  transmitted  are  seven  track,  800  BPI,  unlabelled 
formatted  tapes  with  a record  length  of  120  and  block  size  of  4800. 

The  measured  data  from  the  eight  flight  test  conditions  studied  are 
written  on  two  tapes;  one  containing  five  transient  input  '.onditions  and 
one  containing  three  frequency  sweep  sinusoidal  input  conditions.  For 
each  flight  condition  the  data  appear  on  tape  as  follows: 

FLIGHT  CONDITION  RECORD  (215,  20A4) 


WORD  1 

WORD  2 

WORDS  3 THROUGH  22 

NSEQ 

NODR 

LBL  (20) 

Each  word  in  the  test  data  records  is  a four-byte  integer. 

NSEQ  - Sequence  number  of  the  flight-test  condition  as  it  appears 
on  the  tape 

NODR  - Number  of  data  records  for  the  flight  test  condition 

LBL  (20)  - Alphanumeric  description  of  the  flight  test  condition  (80 
characters) 

TEST  DATA  RECORDS  (REPEATED  NODR  TIMES)  (10E12.6) 


WORD  1 

WORD  2 THROUGH  58 

TIME 

RESP  (57) 

Each  word  in  the  test  data  records  is  a four-byte  real  (floating  point) 
variable . 

TIME  - Time  in  seconds  measured  from  the  beginning  of  the  test 

RESP  (57)  - Recorded  measurements  for  the  57  responses  shown  on  page  453. 

The  f.’i.ht  condition  record  and  corresponding  test  data  records  will  be 
repeated  for  each  of  the  flight  conditions  included  on  the  tape. 
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3.0  Data  Dropouts 

The  recorded  flight  test  data  are  relatively  free  of  data  dropouts.  How- 
ever, the  data  from  two  of  the  flight  conditions,  44.1.5  and  44.3.17, 
contain  dropouts.  The  locations  and  definitions  of  the  dropouts  are: 

Condition  44.1.5  Response  No.  40  - lateral  acceleration  inboard  nacelle 
Data  points  2710  through  2720 

Time  (seconds  fiom  start  of  data)  54.18  through  54.38 
Value  1.0  x lO2^ 

Condition  44.3.17  Response  Nos.  55,  56  and  57  - a , $ and  velocity 
Data  points  548  through  583 

Time  (seconds  from  start  of  data)  10.94  through  11.64 
Value  - 1.0055  x 1042 


TABLE  29.  RECORDED  FLIGHT  TEST  MEASUREMENTS 


! NO. 

ITEM 

LOCATION 

! 

! 1 

Vertical  Bending  Moment 

BS760 

1 2 

Vertical  Bending  Moment 

BS12.22 

! 3 

Vertical  Bending  Moment 

BS1412 

I 4 

Lateral  Bending  Moment 

BS760 

■ 5 

Lateral  Bending  Moment 

BS1222 

6 

Lateral  Bending  Moment 

BS1412 

7 

Vertical  Bending  Moment 

LWS222 

8 

Vertical  Bending  Moment 

LWS820 

9 

Vertical  Bending  Moment 

LWS974 

10 

Chordwise  Bending  Moment 

LWS222 

11 

Chordwica  Bending  Moment 

LWS820 

12 

Chordwise  Bending  Moment 

LWS974 

13 

Vertical  Bending  Moment 

RWS222 

14 

Vertical  Bending  Moment 

RWS820 

15 

Vertical  Bending  Moment 

RWS974 

16 

Chordwise  Bending  Moment 

RWS222 

17 

Chordwise  Bending  Moment 

RWS820 

18 

Chordwise  Bending  Moment 

RWS974 

19 

Vertical  Bending  Moment 

LHTBL56 

20 

Vertical  Bending  Moment 

RHTBL56 

! 21 

Vertical  Bending  Moment 

FS135 

22 

Pitch  Rate 

BS860 

23 

Roll  Rate 

BS860 

24 

Yaw  Rate 

BS860 

25 

Vertical  Acceleration 

BS172 

26 

Vertical  Acceleration 

BS860 

27 

Vertical  Acceleration 

BS1655 

28 

Lateral  Acceleration 

BS172 

29 

Lateral  Acceleration 

BS860 

30 

Lateral  Acceleration 

BS1655 

31 

Vertical  Acceleration 

LWBL565 

32 

Vertical  Acceleration 

LBWL925 

33 

Vertical  Acceleration 

LWS1359 

1 34 

Vertical  Acceleration 

Left  Ext.  Tank,  Nose 

35 

Vertical  Acceleration 

RWBL565 

36 

Vertical  Acceleration 

RWBL925 

37 

Vertical  Acceleration 

RWS1359 

38 

Vertical  Acceleration 

Right  Ext.  Tank,  Nose 

39 

Vertical  Acceleration 

Left  Inboard  Nacelle 

j 40 

Lateral  Acceleration 

Left  Inboard  Nacelle 

! 41 

Vertical  Acceleration 

Left  Outboard  Nacelle 

! 42 

Lateral  Acceleration 

Left  Outboard  Nacelle 

i 43 

Angular  Acceleration,  Y-Axis 

LWS540 

44 

Angular  Acceleration,  Y-Axis 

LWS1359 

45 

Angular  Acceleration,  Y-Axis 

RWS1359 

46 

Delta  Inboard  Aileron 

- 

47 

Delta  Rudder 

- 

453 


UNITS 


In. -LBS 


In. -LBS 
Deg/Sec 
Deg/Sec 
Deg/Sec 


® 0 
Rad/Sec^ 

Rad/Sec^ 

Rad/Sec^ 

Deg 

Deg 


TABLE  29.  RECORDED  FLIGHT  TEST  MEASUREMENTS  (Concluded) 


LOCATION 


Delta  Elevator 
Delta  Vertical  Canard 
Delta  Left  Horizontal  Canard 
Delta  Right  Horizontal  Canard 
Delta  Outboard  Aileron 
Delta  Outboard  Flaperon 
Delta  Inboard  Flaperon 
Angle  of  Attack 
Angle  of  Side  Slip 
Forward  Velocity 


Gust  Boom 
Gust  Boom 
Gust  Boom 


MEASUREMENT  SIGN  CONVENTION 


Vertical  Bending  Moments 


Lateral  Bending  Moment 

Chordwise  Bending  Moment 
Pitch  Rate 
Roll  Rate 
Yaw  Rate 

Vertical  Acceleration 
Lateral  Acceleration 
Angular  Acceleration,  Y-Axis 
Control  Surface  Deflection 
Angle  of  Attack 
Angle  of  Side  Slip 


Fuselage,  + Tail  Up 
Wing,  + Tip  Up 
Horizontal  Tail,  + Tip  Up 
Fuselage,  + Tail  Right 
Fin,  + Tip  Right 
Wing,  + T.  \ft 
+ Nose  Up 
+ Left  Wing  Up 
+ Nose  Right 
+ Up 

+ Right* 

+ Leading  Edge  Up 
+ Trailing  Edge  Down  or  Left 
+ Nose  Up 
+ Nos<s  Lsft 


*Lateral  acceleration  at  BS  1655  is  left 
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